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The bulk viscosity of three-flavor color-superconducting quark matter originating from the non- 
leptonic process it-|-s^it-|-ciis computed. It is assumed that up and down quarks form Cooper 
pairs while the strange quark remains unpaired (2SC phase). A general derivation of the rate of 
strangeness production is presented, involving contributions from a multitude of different subpro- 
cesses, including subprocesses that involve different numbers of gapped quarks as well as creation 
and annihilation of particles in the condensate. The rate is then used to compute the bulk viscosity 
as a function of the temperature, for an external oscillation frequency typical of a compact star 
r-mode. We find that, for temperatures far below the critical temperature Tc for 2SC pairing, the 
bulk viscosity of color-superconducting quark matter is suppressed relative to that of unpaired quark 
matter, but for T > 10~^ Tc the color-superconducting quark matter has a higher bulk viscosity. 
This is potentially relevant for the suppression of r-mode instabilities early in the life of a compact 
star. 
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I. INTRODUCTION 

m I 

^ \ Compact stars consist of extremely dense matter, part or all of which may be deconfined quark matter [ij. Quark 
On , matter that is sufficiently dense and cold is expected to be in a color-superconducting state , characterized by Cooper 
pairing between the quarks, mediated by their attractive strong interaction. This is analogous to the formation of 
electron Cooper pairs in a superconducting metal However, unlike electrons, quarks occur in many varieties: 
' there are three colors and in a compact star core we expect to find three flavors, so a multitude of pairing (and thus 
\^ , symmetry breaking) patterns and hence different phases is theoretically possible (3] • It is an open problem of current 
' research to determine the ground state out of these many possibilities for quark matter inside a compact star. 

One approach to this problem is purely theoretical, using perturbative QCD as well as more phenomenological 
^ . models in order to compare different possibilities and find the one with the lowest free energy. This approach yields 
' a definite answer for very large densities (more precisely for quark chemical potentials ^ much larger than the masses 
^ , of the up, down, and strange quarks). In this situation, quark matter is in the color-flavor locked (CFL) state [H]- 
However, in nature, i.e., in the interior of a compact star, the quark chemical potential is expected to be of the order 
of 500 MeV or less, which requires the mass of the strange quark rus to be taken into account. Moreover, further 
constraints on the matter are imposed by electric and color charge neutrality as well as by equilibrium with respect to 
the weak interaction. All these conditions impose a stress on the pairing, no matter which pairing pattern is chosen 
0. Therefore, besides the CFL phase, many unconventional color-superconducting phases have been proposed, e.g., 
d ' phases that break rotational Q or translational [1] symmetries. 

The second approach to determining the ground state of dense quark matter in a compact star is the study of possible 
astrophysical effects of color superconductors. In particular, the ultimate goal is to compare observational data with 
theoretical predictions for color superconductors in order to exclude or confirm proposed phases. Moreover, one would 
like to use observations to distinguish between quark matter and more conventional matter inside the star, such as 
nuclear or hyperonic matter. Promising physical properties that are distinct for different color-superconducting phases 
are related to transport properties such as thermal conductivity , neutrino emissivity [l^, [isl . [l^ . [isj , magnetic 
fields |T3| , and viscosity. 

In this paper, we calculate the bulk viscosity of color-superconducting quark matter in the 2SC phase. (For a 
recent study of the bulk viscosity in spin-one color superconductors, see Ref. [l^; for hybrid/CFL matter see Ref. 
[l7| . ). In 2SC quark matter, up and down quarks form Cooper pairs with each other while the strange quarks 
remain unpaired. Although constraints of charge neutrality tend to disfavor the 2SC phase [isl . [Toj . it is found in 
phase diagrams calculated from NJL models, for example in Ref. [l^l (see their Figs. 2 and 6) and in Ref. [2l| (see 
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their Fig. 7). In our calculations wc shall neglect the effect of electric charge neutrality, so we treat the up and down 
quarks as having the same chemical potential and hence the same Fermi momentum. This has a quantitative effect, 
which we expect to be of order ml/ fx^ , on the available phase space for the quark momenta. We therefore expect our 
results to deviate only slightly from the results for electrically neutral 2SC matter. 

Bulk viscosity, as well as shear viscosity [2^ . is relevant to vibrational and rotational modes of a compact star, 
whose frequencies are typically in the kHz range [2^. In particular, for a rapidly rota ting star there is the possibility 
of rotational "r-modes" that arc unstable with respect to gravitational radiation [23, HH, [H, [13] ■ This would lead to 
a drastic decrease of the rotation frequency of the star, which is incompatible with observations of rapidly rotating 
|)i^sars with periods of the order of nis. One mechanism to damp the instabilities is a nonzero viscosity of the fluid 

In general, bulk viscosity characterizes the response of the system to an externally imposed oscillating volume 
change. If the volume compression (and expansion) induces a change in the chemical composition of the matter, i.e., 
if it drives the system out of chemical equilibrium, then microscopic processes will try to re-equilibrate the system. If 
these processes occur on a timescale comparable to that of the driving oscillation then power dissipation will occur. 
Because the strange quark is massive, a volume oscillation causes a difference between the chemical potential of 
the s quark and that of the u and d quarks. The relevant microscopic rc-equilibration process is the weak process 
w + (i ^ u + s [2^ which produces (annihilates) strangeness. (Leptonic processes are suppressed by their smaller 
phase space. However, we shall discuss the possibility that they might still yield a significant contribution to the bulk 
viscosity.) Strong processes can be neglected since they occur on much shorter time scales than the scale given by the 
oscillation of the star. This is in contrast to the context of heavy ion collisions where gluonic processes are important 
for the bulk viscosity [s^j- Moreover, we do not consider a potential contribution from the Goldstone boson associated 
with the spontaneous breaking of U{\)a because this symmetry is explicitly broken at densities relevant for compact 
stars. Other (approximate) global symmetries arc not broken in the 2SC phase, hence there arc no other (pseudo) 
Goldstone bosons. 

In order to compute the rate of the process u + d ^ u -\- s one has to consider the corresponding collision integral. 
In the case of a superconductor/supcrfluid this integral is considerably more complicated than in the unpaired phase. 
Some of the complication originates from the broken number conservation, allowing not only processes that consist of 
two ingoing and outgoing particles but also processes with three (four) particles coalescing into one (zero) particles 
and one (zero) particles decaying into three (four) particles. In these processes, particles are created from or deposited 
into the condensate. This effect is also present in more conventional superfluids such as ^He [3lll32|. Starting from the 
kinetic equation, we present a detailed derivation of the collision integral which contains all these possible processes. 

An additional complication arises from the nature of the 2SC phase. The process ui+d ^ U2 + s actually contains 
several subprocesses according to whether the incoming and outgoing states are gapped (up or down quarks that 
are red or green) or ungapped (all the other species). All these subprocesses are included naturally in our formal 
derivation. We encounter processes where zero, one, two, or three of the quarks ui, d, U2, s are gapped. Therefore, the 
2SC phase is an interesting system to study the more general question of the quantitative effects of the energy gap for 
the rate of the microscopic process. We compute the rates for arbitrary temperatures T (up to the superconducting 
transition temperature), i.e., wc go beyond the low-temperature regime, where the effect of the energy gap A is an 
exponential suppression of the rate. 

A final subtlety is that the bulk viscosity is not a monotonic function of the rates of the equilibrating processes. 
The bulk viscosity is highest when those rates are closest to the time scale of the oscillation that is being damped. 
As we will see, this means that, depending what frequencies are physically relevant, a phase in which equilibration 
proceeds more slowly may have a higher bulk viscosity than one in which it proceeds quickly. Thus the 2SC phase, 
whose equilibration rate is always slower than that of unpaired quark matter, will turn out to have a larger bulk 
viscosity at kHz frequencies over a wide range of temperatures. 

The paper is organized as follows. In Sec. |TT] wc define the bulk viscosity and point out an analogy to the power 
dissipation in an electric circuit with a capacitance and a resistance. The derivation of the rate of the process 
u + d^u + s in 2SC quark matter is presented in Sec. IIIII At the end of this section, in Sec. IIII El we discuss the 
physical content of the derived expression. In Sec. IIVI the rate is evaluated in full generality and several limit cases 
(unpaired phase, zero-temperature limit) as well as several subprocesses and their contribution to the full result is 
discussed. We use these results to compute the bulk viscosity in Sec. fVland conclude with Sec. IVII 

II. DEFINITION OF BULK VISCOSITY 

In this section we derive the expression for the bulk viscosity with respect to the process u+s ^ u+d. This derivation 
is based on respective derivations of unpaired quark matter and variations of it for different microscopic processes 
[33I [H, [H, [3^ . We also point out an instructive analogy to an electric circuit with resistance and capacitance. 
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Before we do so, we notice that the following derivation includes several assumptions about other processes that 
potentially contribute to the bulk viscosity. First, there may be a contribution from the leptonic processes M+e <-> d+v 
and u + e ^ s + I'. In Appendix |A] we present a derivation of the bulk viscosity including these processes in addition 
to the nonleptonic process u + s ^ u + d. The result (|A13p is complicated compared to the result (PT|) . derived in this 
section. Because of phase space restrictions, the leptonic processes can be expected to be slower than the non-leptonic 
one. However, as we discuss in Appendix El their contribution can still be significant, depending on the external 
oscillation frequency. In this paper, we do not compute the rates of the leptonic processes microscopically. Therefore, 
we leave the study of their quantitative contribution to the bulk viscosity for future work. 

Second, the derivation of the bulk viscosity below treats the effect of the process u + s ^ u + d solely through 
the flavor chemical potentials fi^, and Hs, disregarding color degrees of freedom. This treatment is correct under 
the assumption of infinitely fast, color-changing strong processes as we show in Appendix [BJ In this appendix we set 
up a more general derivation of the bulk viscosity showing that in the absence of the strong interaction (or for not 
infinitely fast strong processes) the expression for the bulk viscosity is more complicated. For this aspect, see also 
discussion of the results of Fig. [7l 

We start from the definition of the bulk viscosity [13, S, [13, [1^ 

where we assume the volume of the system to oscillate with angular frequency uj and amplitude SVq about its 
equilibrium value Vb, 

y (t) = Vo + 5V{t) , SV{t) = SVa cos ut . (2) 
The average dissipated power per volume in one oscillation period r = 2tt/u! is 

where p{t) is the pressure. We shall express the dissipated power in terms of 6V{t) and the resulting oscillation in 
6fi(t). Here, Sfi is the difference between the s quark and the d quark chemical potential Sfi = ^ y^d which vanishes 
in chemical equilibrium. To this end. we write the pressure for small oscillations as 

- + (^) , + {t) , + • 

Here and in the following, the subscript denotes equilibrium, i.e., dji = 5V = 0, and Us and Ud are the s- and 
d-quark number densities, respectively. We employ the thermodynamic relation 

and use 



^0 1^ = "dTT- + ^uTj— , (6a) 



-Vo - (n.P) . (6b) 



dV / g \ dn. 

Here we have taken into account that, in the 2SC phase, u and d quarks form Cooper pairs with each other. Therefore, 
Hu and /id depend on both Uu and Ud- In contrast, solely depends on Ug because s quarks are unpaired. See Table U 
for the explicit form of the densities. Next, we express the change in densities in terms of the rate of the microscopic 
process, to be computed in Sees. Illll and llVi 



Sudit) = -Srisit) = / Td[6^iit')] dt' . (7) 
Jo 

Here, Td[6fi{t)] is the number of d quarks per volume and time produced in the reaction u + s u + d. Since both 
directions of the reaction are taken into account, the rate is zero for Sfi = 0. In the present situation, however, a nonzero 
6/1 is induced through the volume oscillation, rendering Td[S^{t)] nonzero. We shall use the linear approximation 

Td[Sti{t)]^XSti{t). (8) 
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Inserting Eqs. ([6|) into (O, the result and ([7]) into Eq. the result into Eq. ([3]), and making use of Eq. ([8]) we find 

where the right-hand side has been obtained via partial integration, we have abbreviated (5u(t) = SV{t)/Vo, and 

B^Up^) +(n//^) ^Up-) . (10) 

Finally, we have to eonneet the change in eheniieal potentials to the volume oscillation and the microscopic rate, 

dS^ f d6^i\ dV f dS^\ dn^ f d5^\ drig 



dt y dV ) Q dt ' \ dud J Q dt ' dug J ^ dt ' ^ ^ 

where the change in densities is due to the microscopic process, dud/dt — —dris/dt — Td- Consequently, 

dSn ^ dSv c- , <. / „x 

^-B—-,SKt), (12) 

where we abbreviated the characteristic (inverse) time scale set by the microscopic process by 

7 = CA (13) 



where 

dud/Q \dn 

The equation for the power dissipation ([9]) in terms of the periodic external quantity 6v{t) and the response of 
the system S^{t) together with the differential equation for 6iJ.{t) can be mapped exactly onto the analogous 
situation in an electric circuit. In this case, the external quantity is an alternating voltage U{t) while the response 
of the system is the current I{t). For a circuit with a resistance R, an inductance L and a capacitance C", the 
corresponding equations are 

U{t) ^ i{t)R + Li{t) + l^. (15b) 

Obviously, Eq. (|15b[) has the same form as Eq. and we can identify "voltage" , "current" , "resistance" , "induc- 
tance" , and "capacitance" in the color superconductor, 

1 B 

U(t)^5v(t), I{t)~^Sn(t), R~>—, L^O, C'^—. (16) 

B 7 

This analogy with an electric circuit is not only of a mathematical nature, we can also use it to understand the 
underlying physics. This is useful because the behavior of the bulk viscosity, unlike the shear viscosity, is somewhat 
counter-intuitive if one applies the daily life picture of a "viscous" fluid. Therefore, let us interpret the analogies in 
Eq. ([in]). First, the "resistance" of the system is proportional to 1/_B with B given in Eq. PT7)1 . In particular, the 
resistance is infinite for vanishing B. This means that there is no response of the system to the external force, just 
as there is no current in a circuit with infinite resistance. The reason is that a vanishing B means that there is no 
difference of d and s flavor in the response to the change in density induced by the volume oscillations. Hence the 
system is still in chemical equilibrium and no energy is dissipated. A finite resistance occurs if there is a difference 
in the dispersion relations of d ans s quarks, such as a different mass or different energy gaps. Only then does the 
system respond with an oscillating Sfi(t) and both the power dissipation and the bulk viscosity are non- vanishing. 

The "capacitance" is proportional to I/7, where 7 is the inverse time scale that characterizes the weak process. A 
slow process has a small 7, hence a large capacitance, meaning that the system can store a large amount of energy. 
This energy is only slowly released by producing d quarks (or s quarks, depending on the sign of Sfi) and thus re- 
approaching chemical equilibrium. This is just like a large capacitance which can store a large amount of electrical 
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charge. On the other hand, a fast process has a large 7 and the capacitor can be discharged quickly (i.e. chemical 
equilibrium can be reached quickly). 

Let us now continue with our derivation. We write the response of the system as 

5fi{t) ^Rc{6fiae"^^) (17) 

with a complex amplitude SfiQ. This expression is inserted into the differential equation (jl2p . We obtain a phase lag 
(j) of (5/i(<) with respect to Sv{t), 

= arctan — . (18) 
u; 

The amplitude of the oscillation of the chemical potential is 

l^.ol = f^^|^, (19) 



and the power dissipation is 



2 n2\,,2 

(20) 



\ dt / 2 Vo 2 V Vb / 72 + ^2 ■ 

With the definition ^ we arrive at the following expression for the bulk viscosity, 

7 

C = 7^^^: (21) 

G 7^ + oj^ 

This expression, together with the definitions ([TU]) . ([13]), and (O, is the main result of this section and shall be used 
in Sec. El 

In Tabic [H we show the explicit values for the coefficients B and C in unpaired and 2SC quark matter at zero 
temperature, derived from the quark densities as functions of the chemical potentials. The densities can be computed 
from the respective free energy in Ref. [l^. In order to obtain the partial derivatives appearing in B and C in the 
2SC phase, one has to invert the Jacobian of the two-dimensional function (nu^rid) ~ {'nu{^J■u^ ^J'd),nd{^J-u, ^J'd))■ We 
have neglected terms in higher orders of and the pairing gap A. Moreover, we expect the difference in chemical 
potentials Siy = {fj-u — IJ-d)/'^ to be of the order of ml/ [l3. Therefore, we have also dropped higher order terms in 5v. 
In the following calculation of the microscopic rate, we shall, for simplicity, consider the limit /i„ ^ fid- It is important 
to note that even in the limit A = the value of the coefficient C is different in the two considered phases. The reason 
is the "locking" of chemical potentials of the pairing fermions. In other words, gapped quasiparticlcs in the 2SC phase 
look as if their chemical potential (in the absence of pairing) was v = [jiu + /^d)/2. This effect, the difference in C for 
unpaired and 2SC matter, is weakened by a nonzero temperature. We have checked numerically, however, that the 
difference remains sizable for all temperatures up to the transition temperature. Thus, for simplicity, we shall use the 
zero-temperature values from Table U as an approximation. 

We have checked that electric neutrality has no sizable effect on the coefficients B and C . Requiring local electric 
neutrality at all times yields a constraint for the densities and chemical potentials via the condition 2n„ — — ns = 0. 
Implementing this condition modifies the derivation for the bulk viscosity. However, to leading order the value of 
B^ jC turns out to be the same as in the non-neutral case. 



III. DERIVATION OF THE RATE FOR THE PROCESS u + s <-> u + d 

In this section we derive the rate r^, see Eq. ([7]), as a function of b[i^ i.e., for a system in chemical non-equilibrium. 
In equilibrium, = and = 0, since in this case the processes u + s u + d and u + d —f u + s have the same 
rate and thus the net rate r^; vanishes. We shall assume Sfi — fig ~ f^d > 0, such that we expect a net production of 
d quarks, Td > 0. The main result of this section is given by Eqs. (|47p . ((50)) . and ((5T|) . 



A. Formalism 

We start from the kinetic equation which is derived from the Kadanoff-Baym equation [s^] after applying a gradient 
expansion. For similar calculations within this formalism see Refs. (isl. [l^. Issl [39t . Close to equilibrium and employing 
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unpaired 


2SC 


riu 




+ 27/3/3 + I.) /tt^ 


rid 




(Md/3 + 2//V3 + A27y)/7r^ 




illl-K^ [1 - 3771^/(2^?)] 


/i^/TT^ [l-3m./(2^7?)] 


B 


m^/(3/id) 


m?/(3Md) + 2AV(9i/) 


C 


2^V(3/il) [1 + m2/(4/i^)] 


TrV/il [1 + m2/(6/il) - AV(18Md) - 25/^/(3^^)] 






mi/(97r2)+4A2m^/(277r2) +4AV(8l7r2) 



TABLE I: Quark densities tiu, ti^, tIs and coefficients B and C for unpaired and 2SC quark matter for massless it and d quarks 
and a strange mass ttis <S iJ,u,Hd,Hs. We liave abbreviated v = {jj.u + Pd)/2 and Sv = {fiu — /id)/2. Tlie energy gap due to 
pairing of it and d quarks is denoted by A. For instance, tlie it quark density in tlie 2SC pliase is given by tfiree contributions: 
unpaired blue quarks yield a term ^/{Stt^); paired red and green quarks yield two terms, one contribution from a Fermi sphere 
filled up to the average chemical potential u plus condensation energy proportional to A^. In the last line we give B^/C to 
lowest order in rris and A, which determines the maximum bulk viscosity. 



the closed-time-path formalism we may use the kinetic equation 

d 



z^Tr[7o5<(Fi)] = -Tr[5>(Pi)S<(Pi) - S>(Pi)5<(Pi)] 



(22) 



Here and in the fohowing, four-momenta are denoted by P = (po,p)- The quantities 5^'^ and E<'> arc the fermion 
propagators and self-energies, respectively. More precisely, they represent the d-quark propagator and the contribution 
to the d-quark self-energy originating from the (left) diagram shown in Fig. [1] Then, the left-hand side of Eq. (|22p . 
integrated over Pi, basically is another way of writing F^, while the right-hand side yields the collision integral for 
the process u + s u + d. Hence, in the following, we shall compute the rate by evaluating the right-hand side of 
Eq. (P^ . To this end, we first define S"^'^ and for the 2SC phase (remainder of this subsection and Sec. IHI B)) 

and then insert these definitions back into the right-hand side of Eq. ([22)1 to obtain the collision integral (Sees. IHI CI 
andHHD]). 

Since we are describing a superconducting system, the fermion propagators are matrices in Nambu-Gorkov space. 



5<'>(P) 



G<'>(P) P<'>(P) 



p_;-"(p) Gi'^(p) 



(23) 



The explicit form of the "normal" and "anomalous" propagators G^'^ and P^'^ shall be given in the next subsection 
for the 2SC phase. The relevant part of the c?-quark self-energy is given by the diagram in Fig. [1] (the right diagram 
in this figure is merely a reminder of our notation of the momenta). 



(Pi) 



(27r)4 



^<^>(^4)rL + H>-<(Pi-P4) 



(24) 



where Mw is the W^-boson mass and F,^^ _|_ is the vertex of the subprocess d ^ u + W . The respective contribution 
to the VF-boson polarization tensors are given by 



n<^>(g) 



(27r)- 



-Tr[Fi:,,+ 5>-<(P2 + Q) F.L._ 5<->(P2)] , 



where F(^^ _|_ is the vertex of the subprocess u + W 



s. In Nambu-Gorkov space, the vertices are 



(25) 



ud/us,:t 



d/us / 7^(1 - 7^)T„rf/„ 



sm Ow 







-7^(l+7'^)^nd/n 



with the entries of the CKM matrix Vud and Vus and the weak mixing angle dw, and 
^0 1 

Tud.+ ^ 

\0 





(0 














') 






(0 





Tud- = 


1 





=) 


, Tus,+ = 











1 'us, 





















^0 





«) 












(26) 



(27) 



are matrices in flavor space (with indices ordered u, d, s). 
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B. Quark propagators in the 2SC phase 

We now specify the normal and anomalous propagators G^'^ and F^'^ . We shall consider the simplest situation 
of massless quarks. In this case, the propagators in the 2SC phase are 

G<>{K) ^jZVrGt->{K)i'ht , F<'>{K) ^ J^hF<:>{K)^''A^ . (28) 

The index r labels the different quasiparticle excitations, which are explained below and shown in Table [TTl The 
quantities = (1 T 7^7 • k)/2, with k = k/fc, are energy projectors, {Ii)jk = {Ji)jk = ^i^ijk < 1,2,3) are 

matrices in flavor and color space, respectively, and 

G^,,(X) = -2Tn fiek,r) 6{ko ± Mr - e^.r) + [1 - /(e^,,)] S{ko ± Mr + efe.r)} , (29a) 



< 

±.r 



(K) = -2^1 [1 - !{ek,r)\ S{ko ± Mr - £k,r) + Bl^. f{ek,r) S{ko ± Mr + £k,r)} , (29b) 

F^r{K) = 27ri -^|/(£fc,r)(5(fcoT Mr - £fc,r) - [l-/(£fe,r)]<5(fco T Mr + £fc,r)|, (29c) 

ZSk.r ^ J 



k.r 

[K) = 27ri^{[l-/(efc,0]<5(A:o±Mr-£fc,r)-/(£fc,r)<5(fco±Mr + £fc,r)}. (29d) 



2efc.r 

Here, we have introduced the Bogoliubov coefficients 



^.^^£^iA^_ (30) 



Moreover, / is the fermion distribution function. For the right-hand side of the kinetic equation, we shall approximate 
/ by its equilibrium form. 

We have neglected the antiparticle contribution to the propagators. For simplicity, we assume Mm = Md = This is 
somewhat unrealistic since it violates overall electric charge neutrality. However, a more realistic scenario Mu Md 
would lead to some additional (notational) complication and does not alter our results qualitatively. With this 
simplification, we have 3 different quasiparticle dispersion relations, Ek.n with r = 1, 2, 3. They arc shown in Table Ull 
with A being the energy gap induced by Cooper pairing of u and d quarks. The condensate picks a certain direction 
in color and flavor space, here J3 and J3, respectively, leaving one color of u and d quarks as well as all s quarks 
unpaired. Formally, each of the three modes lives in a subspace of the 9-dimensional color-flavor space, which is given 
by the corresponding projectors Vi occurring in the normal propagators G'^'^{K). The projectors are shown in the 
last column of Table[TTl The dimensions of the subspaces are 4, 2, 3, respectively. In the anomalous propagators, only 
the gapped modes, r — 1, contribute. 
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r 


quarks 


dispersion 


projector 


1 


ru, rd, gu, gd 




Pi = lUi 


2 


bu, bd 


£k.2 = \k - ^\ 


V2 = ilii-Ji) 


3 


rs, gs, bs 


£k/i = fc — ^is 


p3 = 1 - li 



TABLE II: Quasiparticle modes in the 2SC phase for ultra-relativistic quarks and jj.u — fJ-d = 



C. W-boson polarization tensor 

Next, we compute the contribution to the M^-boson polarization tensor that occurs in the rf-quark self-energy, Eq. 
24|) . From its definition ([25]) we obtain after inserting Eq. ([23| and performing the trace over Nambu-Gorkov space 



n<^>(Q) = / ^{W(1-7')wG>'<(^2 + Q)7"'(1-7')t„..-G<'>(P2)] 

o sm (7vi/ J l^""; ^ 

+ Tr[7'^ (l + 7^) r„,._ G>'< (P2 + 0) 7^(1 + 7') rus,+ G< > (P2)] 

^ Trini + 7^) r„,,_ F>'< (P2 + Q) 7''(1 - 7'') P<'>(P2)] } . (32) 

We assume without explicit proof that the second (fourth) term yields the same contribution as the first (third) term 
because this is just the charge-conjugate counterpart. Next we use Eq. (|28p and perform the trace over color and 
flavor space. We make use of 

(33a) 
(33b) 
(33c) 
(33d) 

Because of Eq. (|33d[) . the anomalous propagators yield no contribution. Then, with P3 = P2 + Q and denoting 

T^^-^ip, k) ^ Tr[7^(l - 75)70Ap7''(l - 7^)7° A^] , (34) 

we obtain 

ny(g)^-2» / ^T'^-(p3,p2)[2G>;3<(P3)G<;>(P2) + G>;<(P3)G<;2>(P2)] . (35) 

Inserting the propagators from Eqs. (|29p into this expression and using 1 — /(e) = /(— e) for the equilibrium function 
yields 





= , 








= 0, 




-(1- 


Tus, + V3 


= '''tis,-!- 




0, 




= 0. 







yw — J (27r)3 

6263 " ^ ' 

X {2-Bp^,3-Sp',i /(±e3£p3,3)/(Te2£p2,i) '^(go - e3£rp3,3 + e2£p^,i + S^i) 

+ Bpl^iBll -^ /(±e3£p3,3)/(Te2£p2,2) ^(^o - e^ej,.^:!, + e2ep.,^2 + (5^)} , (36) 

where the sum runs over 62, 63 = ±1, and the upper (lower) sign corresponds to 11^ (n^). Remember that the quark 
loop in the polarization tensor contains a u and an s quark. Therefore, both contributions in curly brackets contain 
the third mode, r = 3, describing an (ungapped) s quark. For the u quark, there are two options. The first term in 
curly brackets describes a gapped u quark, r = 1. The factor 2 in front of this term originates from color degeneracy: 
Both red and green u quarks are gapped. The second term in curly brackets describes an ungapped blue u quark. 
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D. Collision integral 

Upon inserting Eq. ((24)) into Eq. (|22|) . the kinetic equation becomes 

-rL,-^>(^4)0,+^<(Pi)n<,(Pi - P4))] . (37) 

Hence, in order to evaluate the right-hand side of this equation, we need the traces 

8sin 0w ^ 

+ Tr[G>(Pi)7^(l +75Kd.+G<(P4)7"(l + 7')T„d,-] 
- Tr[P>(Pi)7^(l +75)r„,.+P<(P4)7'^(l - 7')r„rf,+] 
-Tr[P>(Pi)7^(l-75K,._P_<(P4)7'^(l + 7'')r„rf,_]}, (38) 



Tr[r^:, _5>(P4)0.+5<(Pi)] = ^!^|Tr[7^(l - 7')r,,,_G>(P4)7''(l - 7')r,,.+G<(Pi)] 

8 sm I- 



and 

i sin^ Ow 

+ Tr[7^(l+75)T„d,+G>(P4)7''(l + 7')^.d,-Gf(Pi)] 

- Tr[7^(l - 75)r„d,_P>(P4)7''(l + 7')T„d,-P< (Pi)] 

- Tr[7^(l +75)r„,.+P>(P4)7''(l - 7')r„rf,+P<(Pi)]} • (39) 

Again, in both expressions, the first (third) and second (fourth) traces yield identical contributions. For the traces 
over color and flavor space we make use of 

T„d,±Pi = Tud,±Jl , (40a) 

Tud,±V2 = T„d,±(l- Jl), (40b) 

Tud^±V3 = 0, (40c) 

Tud,+hTud,+h = If - I (40d) 

Tud^-hTud,-h = ll-l- (40c) 

This shows that, in contrast to the W^-boson polarization tensor, the traces containing anomalous propagators do not 
vanish. We obtain 

Tr[5>(P:)E<(Pi)] = I ^{ [2G^,i(Pi)G< 4(P4) + G> 2(Pi)G< 2(P4)] r'^'^(P4, Pi) 

4M^ sm Ow J (^i") 

+ 2P>i(Pi)P<i(P4)Zi'^'^(pi,P4)} n>jPi - P4) , (41) 

and 

Tr[E>(Pi)5<(Pi)] = ^ ^J^^hn I ^{ [2G> i(P4)G< ^(Pi) + G> 2(P4)G< ^(Pi)] r'^''(p4, Pi) 

+ 2P>i(P4)P<i(Pi)Z^^''(p4,pi)} n<,(Pi - P4) , (42) 
where we defined the trace over Dirac space, 

Z^'^"^(p, k) s Tr[7^(l - 7^)75 A+7"^(l + 7^)7' A^] • (43) 



Again, the different contributions in curly brackets in Eqs. (|4T|) and (|42)) are easy to understand. Remember that 
these expressions describe the subproccss d W~ + u. Therefore, one contribution comes from the gapped u and 
d quarks, r = 1. This mode also yields a contribution from the anomalous propagators. It describes red and green 
quarks, hence the prefactor 2. The second contribution comes from the ungappcd mode, r = 2, which describes blue 
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quarks and which, of course, yields no anomalous contribution. Due to color conservation of the weak interaction, no 
mixing of the two modes is possible. 

Next, we perform the integration over p^Q. Moreover, we shall integrate the kinetic equation over d-quark energies 
Pio and momenta pi, because we are interested in the total rate, not in the rate as a function of momentum. Therefore, 
we need the results (derived by inserting the propagators from Eqs. ((29)) into Eqs. ((4T|) and p2)) ) 



J-oo 27r ^ ^ 4M4^sin2^^y Z_.7 (27r)3] 



b;i,b;i, t'-(p4, pi) + w^'iPuM 



X /(eiepi,i)/(-e4ep^.i)n> (eiEpi,! - e4ep^s,Pi - Pi) 



-^pL2^pi,2'^^''(P4,Pi) /(eiepi,2)/(-e4£p4,2)n>^(ei£p,,2 - e4ep,,2,Pi - P4) L (44) 



and 



r %i£Tr[i]>(P05<(Pi)] = -^^^^y: I 



^pLl-ga,l^^^(P4,Pl)+ Z^^-(P4,P1) 



X /(-eiepi4)/(e4ep^4)n<^(eiepi_i - e4ep^_i,pi - P4) 

+ 5pJ,25pi,2'^^''(P4,Pi) /(-ei£pi,2)/(e4£p4,2)n<^(ei£p,,2 - e4£p,,2,Pi - P4)| • (45) 

Then, in order to obtain the right-hand side of Eq. (|22|) (integrated over Pi), we have to insert Eq. ([36]) into Eqs. 
(UU and (pS)) and the result into Eq. ([221) • One needs the explicit results for the Dirac traces and contractions of 
the form T^^T^^ and U^'^T^„. They are given in Appendix [Cl As shown in this appendix, the contractions of the 
form Z^/'^'^T^iy, i.e., the ones originating from the contribution of anomalous propagators, are complex. However, after 
the angular integration in the collision integral, the imaginary parts vanish, see Appendix ID 21 Therefore, we shall 
continue with only the real parts of these contractions. 

As mentioned above, the left-hand side of the kinetic equation ((22l) stands for Td, the change of the total number 
of d quarks per time and volume (as given by the right-hand side). In order to define properly and extract the 
correct numerical factor, we integrate the left-hand side over Pi and project onto the d-flavor subspace. The latter is 
done by inserting the projector Qd = diag(0, 1,0). Moreover, the change in density is obtained upon multiplying by 
the matrix T3 = diag(l,— 1) in Nambu-Gorkov space (see for instance Ref. [s^l)- Then, using the propagators from 
Eqs. ([29ll . we obtain 

'it I 0Tr[7oQ.r35<(Pr)] = -4^ | ^ {2B;^,/(£p,,) + 25+ [1 - /(Sp^i)] 

+ Sp,,2/(£w.2)+S+,2[l-./(£w,2)]} 

where we have used ^ -I- B^.^ ^ = 1. Eq. (|46|) is the microscopic definition of Td which has been used implicitly 
in Sec. |TT1 The factor 4 arises after performing the trace; it originates from the Nambu-Gorkov doubling and spin 
degeneracy. 

Putting the results for the left-hand and right-hand sides together and dividing both sides by 4 yields 

Td = 4:Tf^^ + 2Tf^^ + 2rf^^ + Tf^^ + 4?^"^ + 2^^231 ^ (-47-) 
We shall explain the physical meaning of the several terms in the next subsection. Their explicit form is 

d^Pi f d^p2 f d^p3 f d3p4 



10Q 4^2ta2t^2 M Pi f P2 M P3 M P4 c/ , x 

i28^G,K.K. E y (^'5(pi+P2-P3-P4) 

61626364 

X (1 - pi • p2)(l - Pa • P4) Bl'Bl'B^^'B2^S{eiei + e^e^ - e^e^ - 6424 + 

X [/(eiei)/(e2e2)/(-e3£3)/(-e4e4) - /(-eiei)/(-e2e2)./(e3e3)/(e4e4)] , (48) 
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and 

Pnrarari ^ RA-,r4r:2„ 1/2,1/2 /" 

X [(1 - P4 • P3)(l - pi • P2) + (1 + pi • P3)(l + P4 • P2) - (1 + Pi • P4)(l + P3 • P2)] 

X - — B^^B^'' - — d(eiei + 6262 - 6363 - e4e4 + Sfi) 

ZSi Z£4 

X [/(eiei)/(e2e2)/(-e3e3)/(-e4e4) - /(-eiei)/(-e2e2)/(e3£3)/(e4£4)] , (49) 

with the Fermi coupHng constant Gp- We have introduced an additional integration over pa and the corresponding 
i5-function for momentum conservation. This was impHcitly present in the previous expressions. Moreover, we have 
abbreviated the quasiparticle energies with momentum pi and excitation branch by e, = Ep^^n and the BogoUubov 
coefficients by B^' = Bp'. ^. . We may now perform the angular integration. This calculation is identical to the unpaired 
phase [iO], see Appendix ID] The results are 

p.,r..3.4 ^ Gj^V^^Vl ^ r i[p,^p^^p^^p^)Bl-B^^^Bl^Bl-5(e,e,+e^e2~e^e^-e4e4 + 5pi) 

X [/(eiei)/(e2e2)/(-e3£3)/(-e4e4) - /(-eiei)/(-e2e2)/(e3e3)/(e4e4)] , (50) 

and 

pn.2.3.4 ^ ^^YMk ^ ( J(p^,p2,P3,P4)^i32=i33=^^^(eiei+ 6262 -e3£3- 6464 + <5m) 

l^'" 616^:^64 •^P1P2P3P4 2£4 

X [/(eiei)/(e2£2)/(-e3e3)/(-e4£4) - /(-ei£i)/(-e2£2)/(e3e3)/(e4£4)] , (51) 
where the functions / and / are given in Eqs. (|D13[) and where we abbreviated 

j /"OO /"OO /"OO 

dpi / dp2 / dpz / dpi . (52) 

P1P2P3P4 



E. Subprocesses and unpaired phase limit 

The collision integral is the sum of six terms, see Eq. (|47p . which are shown diagrammatically in Fig. [2] (in the 
same order as they appear in the equation). The superscripts at Yd and Td in Eq. (|47p correspond to the momenta 
in the order pi, p2, P3, pa, i-c, they describe the particles d, ui, s, U2, respectively, see diagrams in Fig. [T] The 
combinatorics of the process is constrained by color conservation for the vertices d <-> U2 + W~ and ui + W~ s. 
Consequently, d and U2 as well as s and ui must have the same color. Hence, from Table|ll]we conclude that d and U2 
must be either both gapped (r = 1) or both ungapped (r = 2). This can be seen in Fig. [21 where we have added the 
gap A at the fermion line whenever the corresponding mode is gapped. It can also be seen in Eq. (|47p . where the pair 
of first and last superscripts is either (1,1) or (2,2). Furthermore, the second and third superscript can be (1,3) or 
(2,3), describing a (red or green) gapped u and a (red or green) s or a (blue) ungapped u and a (blue) s, respectively. 
The s quark is ungapped for all colors. Combining these pairs and counting color degeneracies yields the six terms 
with respective prcfactors in Eq. (|47)) and Fig. [2l In the cases where both d and U2 are gapped, there is an additional 
contribution from the anomalous propagators. This is denoted by in Eq. (|47[) and by a condensate insertion in 
Fig. [21 From the diagrams it is clear that there can be no condensate insertion for ui. This would lead to a violation 
of electric charge conservation at the vertex ui + ^ s. Summarizing the result of the combinatorics, we have 
found that the process d + ui <-> s + U2 in the 2SC phase can be considered as composed of several subprocesses, 
where zero, one, two, or three of the participating particles have a gap in their excitation spectrum. 

From Eqs. (jSOj) and (|5T|) we see that the process d + m ^ s + U2 not only contains several subprocesses due to the 
2SC pairing pattern, also every single of these subprocesses is written as a sum over 16 terms, expressed as a sum over 
ei, . . . ,64. To understand the physical meaning of this sum, it is instructive to rewrite it and recover the unpaired 
phase limit. The change of variables described in the following is illustrated in Fig. [3l We may rewrite every single 
sum over ei by observing that for any function F we have 



^ / dp^ B^'F{e,e,) = ^ / dp, Bf-F(-e,£,) , 

et -^0 6, -^0 



(53) 
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FIG. 2: Contributions to the process d + u ^ u + s. The four diagrams in the first row contain 3, 2, 1, and gapped fermions, 
respectively (a gapped fermion is marked by A at the respective lines) . The combinations of quarks of the first two diagrams in 
the first row allow for nonvanishing contributions of anomalous propagators. These are shown in the second row (a condensate 
insertion is represented by a hatched circle). 



0.5 



-0.5 




FIG. 3: (Color online) Thick (red) lines: gapped modes. Thin (blue) lines: ungapped modes. Solid lines: + e. Dashed lines: 
— e. The quantities are given in arbitrary units, /x = 1. This figure illustrates the new variables introduced in Eq. (|54p. which 
are advantageous for the discussion of the various subprocesses, as explained below Eq. (|57|l . 



where the new Bogoliubov coefficients and quasiparticle energies are defined as 



(54) 



Consequently, in Eqs. (|50p and (|5ip we can replace by i?^^' and eiEi by —CiSi. Then, we obtain for the contribution 
from the normal propagators. 



£1626364 PlP^PSPi 

X [f {-eiei)f {-6262)1 (esis)/ {6464) - f {eiei)f {€262)/ {-6363) f {-€464)] , 



(55) 



and for the contribution of the anomalous propagators 

T/2 T/2 
'^F ^ud ^us 



^ d 



E 



^ 264 



,p 5 j -^bl, P2,P3,P4)y^ 52^53^ 7^ (5(eiei+ 6262 -6353- 6464 -(5/i) 

61626364 PiP^PSPi ^^^^ 



X [f {-eiei)f{~e2e2)f (6363) f (6464) - f{eiei)f (6262) f (-6363) f (-6464)] ■ (56) 
Because of = for A = 0, only one of the 16 terms, namely ei = . . . = 64 = 1, survives in the unpaired phase 
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limit and wc immediately obtain for this case 

punp ^ ^^pVudVus / /(P1,P2,P3,P4)<5(P1+P2-P3-P4) 
X {[1 - /(pi - /i)][l - f{p, - ^)]f{p, - IJ,)f{pi f{p, - ^)f{p2 - /i)[l - f{p3 - /is)][l - /(P4 - /i)]} , (57) 

in accordance with Rcf. [4lj. Note that in the representation of Eq. ([5D1) . all 16 subprocesscs contribute even in the 
unpaired phase (of course yielding the same result (j57p ). It is easy to understand that the quantities in Eq. (j54p are 
more natural to describe the unpaired phase because Si = Pi —■ describes particle excitations for all momenta pi. 
In contrast, = \pi — describes hole excitations for pi < fi and particle excitations for pi > ^i. This distinction 
between particles and holes is not possible in the quasiparticle description of the superconducting phase. However, 
also in the 2SC phase the new representation is instructive since it reveals the processes that are only possible in a 
superconductor/superfluid, namely the 15 subprocesses that contain at least one B~ . They correspond to processes 
in which three particles coalesce to one particle (or one particle decays into three) and in which four particles are 
annihilated (or created). This is possible because particle number is not conserved. Particles can be created from or 
deposited into the condensate. Only when there is no condensate, the "conventional" process of two incoming and 
two outgoing particles is the only possible one. In the next section, we shall study the quantitative contribution of 
these "exotic" processes to the total rate. 

As a result of this subsection, let us rewrite the collision integral (|47p using the following more illustrative notation, 

Td = AVl'''-'^' (A, A, A) + 2r;^i'=^ (A, 0, A) + 2VJ (0, A, 0) + Td{Q, 0, 0) + AVl'"'-'''' (A, A, A) + 2?;^^*=^ (A, 0, A) , (58) 

where summation over the respective = ± is implied. Wc shall also use the shorthand notation rci(A,A,A) = 
r^^'^^'^'' (A, A, A) etc. where no confusion is possible. This notation shows explicitly how many gapped modes partic- 
ipate in the respective subprocesscs. The definitions of the new quantities is clear from comparison with Eq. (j47p . 
the sum over ei, 62, 64 originating from Eqs. (|55p and (j56p . Since the s quark (corresponding to 63) is ungapped, 
only terms with 63 = +1 survive, and the superscripts ei, 62, 64 as well as the three arguments of the functions 
Td correspond to the quark flavors d, ui, and U2. It is clear from the definition that the number of gapped modes 
determines how many of the 16 subprocesscs arc nonvanishing. This number can be directly read off from the notation 
in Eq. (|58p . i.e., all e,;'s that are shown explicitly have to be summed over ±1, while the missing Ci are set to +1, 
since the contribution for —1 vanishes. 



IV. EVALUATION OF THE RATE 



The general evaluation of the collision integral has to be done numerically. As a consistency check, let us derive the 
analytical result in the unpaired phase for zero temperature. In this case, all Fermi distributions in Eq. (|57p become 
step functions, 

T/2 T/2 c 

/ ^(Pl,P2,P3,Pl +P2 -Ps) 0(Pl +P2 -Ps) 

^ PxPiP?, 

X [e(pi - /i)6(p2 - f^)Q{fJ-s - P3)9(yu - Pi - P2 + Pa) - 0(m - Pi)©(a* - P2)0(P3 - Ms)0(pi +P2-P3- m)] • (59) 

Since we assumed that Sp, ~ fis — /-i > 0, the second term in square brackets vanishes. (The first three step functions 
in this term imply pi + p2 ~ P3 < 2/i — fig which, for /is > fi, contradicts the fourth one, implying pi + P2 ^ P3 > P-) 
Physically, this means that at zero temperature the process only happens in one direction. This direction is given by 
the sign of Sfi. The first term in square brackets can be evaluated upon using the definition of the function / in Eq. 
(|D13ap . One finds 

g /"Ms rP3 /•P3-P2+/1 

C = ^G|Fi,V^,y dp^J dnj <ip,(pi+P2-P3)'[(Pi+»)'+3Pa(Pi+») + 6pH 

- ^.OlVlVi, („=*.3 ^ _ ^ _L,.,y + _ , ,60, 

which is in agreement with Ref. [4l| . 

We treat all other cases numerically. For the following results we use ji = 500 McV and J/i = 10 McV. We 
assume a transition temperature of Tc = 30 MeV. Via the relation Tc/Aq = /ir, where 7 ~ 0.577 is the Euler- 
Mascheroni constant, this corresponds to a zero-temperature gap Ag ~ 52.6 MeV. We shall employ the following 
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model temperature dependence of the gap, 



A(r) = AoWl 



(61) 



We have to compute the colHsion integrals r(j(Ai, A2, A4) for all combinations of Ai, A2, A4 shown in Eq. ((58)) . Let 
us analyse the various subprocesses before we add them up to the total rate. To this end we pick the first term on 
the right-hand side of Eq. ([55]) . which actually is a sum of 8 subrates. We show the rates of all 8 subprocesses of 
r(i(A, A, A) in Fig. 21 The figure shows that, due to symmetries of the collision integral, some of the resulting curves 
coincide. One obtains four groups, characterized by zero, one, two, and three factors B~ in the collision integral 
(from top to bottom in the figure). In other words, the three lower curves in the right panel describe processes of 
coalescence or decay of particles, while the upper curve corresponds to a "normal" interaction process of two incoming 
and two outgoing particles. The left panel shows that this "normal" process is the only one that does not vanish 
at the transition temperature. Of course, all other processes can only contribute to the rate in the presence of a 
condensate. Although smaller in magnitude, these processes cannot be completely neglected for any temperature 
range since their contribution increases quickly for temperatures decreasing down from Tc- A quantitative analysis 
shows that the seven processes where at least one of ei, 62, 64 equals —1 add up to about 30% of the total rate at 
T/Tc ~ 0.9, and their relative contribution increases monotonically with decreasing temperature. 





FIG. 4: (Color online) Subprocesses r^^*^^*^* (A, A, A) on a linear scale (left panel) and a logarithmic scale (right panel). The 
uppermost (black) line corresponds to the process r^^+ with two incoming and two outgoing particles. All other processes 
include creation or annihilation of particles from the condensate. These processes are only possible in the superconducting 
phase and thus vanish at the critical temperature. 



An evaluation of the 8 subprocesses from Eq. rc;(A, A, A), shows that their contribution is of the or- 

der of the processes (A, A, A). Or, to give another quantitative estimate of these processes, we find that 

Seiesej ^^^''^''''(A, A, A) < F^^^(A,A,A). Therefore, one can neglect these contributions, i.e., we may use 

Frf(Ai, A2, A4) for the almost identical rate Fd(Ai, A2, A4) + rrf(Ai, A2, A4). 

For the rates F^^''^''(A, 0, A) and F^^(0,A,0) there are 4 and 2 subprocesses, respectively. We do not show the 
results for these subprocesses explicitly, since the corresponding plots are very similar to the ones in Fig. 3) 

After having discussed the subprocesses for a given number of gapped fermions, we shall now compare the rates for 
processes with different numbers of gapped fermions. I.e., so far we have picked the first of the terms on the right-hand 
side of Eq. (|58p and analyzed the subprocesses hidden in this term. Now we shall compare the different terms in 
this equation with each other. This is done in Fig. (5] Every curve in this figure corresponds to the rate of a given 
number of participating gapped fermions, where the number of gapped fermions increases from the top to the bottom 
curve. One sees that all rates containing at least one gapped fermion are strongly suppressed for small temperatures. 
This suppression can easily be extracted analytically, see Appendix [E) For < Aq, the leading behavior for small 
temperatures is 

rrf(A,A,A)^e-(2A-*'')/^, Frf(A,0,A)^e-^/^, F,(0, A, 0) ^ g-^^-*'')/^ . (62) 

In Appendix [El also the behavior for > Aq, is derived (although this regime is not relevant in the present context). 
From the suppression due to the gap (j62p we conclude that for small temperatures, we may approximate the total 
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2SC rate by the contribution of ungapped fermions, Td ~ rrf(0, 0,0) In Fig. [6] we see that this approximation 

starts to break down for temperatures T > lOMeV ~ 0.3 T^. 




0.2 0.4 0.6 0.8 1 

T 



Tc 

FIG. 5: (Color online) Contributions to the 2SC rate, separated by the number of participating gapped fermions. The more 
fermions are gapped, the lower is the rate of the corresponding subprocess. The lowest (green) curve is the sum of the 
subprocesses shown in Fig. H) The sum of all lines in this figure is the total rate Fd in the 2SC phase. 




2SC 
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T [MeV] 

FIG. 6: (Color online) Total 2SC rate (solid, blue) obtained from Eq. (|58p compared to the unpaired phase rate FJJ"'' = 
9Fd(0, 0, 0) (dashed, black) and the rate Ftj(0, 0, 0) (dashed-dotted, red) on a logarithmic temperature scale (remember Tc — 30 
MeV). For small temperatures, the result approaches the zero-temperature result (which, in the case of the unpaired phase, is 
given by Eq. (|60|l 'l. The 2SC rate is well approximated by Fd (0,0,0) up to T ~ 0.3 Tc. Above this temperature, effects of the 
gapped modes are not negligible. For all temperatures, the rate of the 2SC phase is smaller than that of the unpaired phase. 
In Sec.|V]we shall see that this is not true for the bulk viscosity. 



V. RESULTS FOR THE BULK VISCOSITY 

In this section, we compute the bulk viscosity of 2SC quark matter as a function of temperature. We first discuss 
the role of the subprocesses in Sec. IV Al before we compare the total bulk viscosity of the 2SC phase with that of 
unpaired quark matter in Sec. IV Bl We make use of the expression for the bulk viscosity derived in Sec. Ull see Eq. 
(PT|) . This equation shows that we need the quantities uj, B, C, and 7 = CA. For the coefficients B and C, see 
definitions pO)) and ((T4|) and Table [H we use a strange mass rris = 100 MeV and a chemical potential of /i = 500 MeV. 
In order to determine 7, we need the results for the rate F^, presented in the previous section. Remember that, in the 
derivation of the bulk viscosity, leading to Eq. (PTjl . we have assumed the rate to be linear in S^, see Eq. This 
linear approximation is, for any given temperature T, valid for a sufficiently small Physically, this means that we 
assume the volume oscillation SVq/Vo to be sufficiently small, because the magnitude of dfj, is proportional to SVo/Vq, 
see Eq. ([T^ . Technically, we have to compute the rate at a sufficiently small 6fi such that we can determine A, 
and thus the characteristic frequency 7, via A = Td/Sfi. (Hence, we cannot necessarily use the numerical values from 
the previous section which have been obtained for a fixed 5fi = 10 MeV.) 
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A. Contribution of subprocesses to the bulk viscosity 

Before we turn to the final results for the 2SC phase, let us discuss the contribution of the subprocesses. As we have 
discussed in Sec. IIVI the rates of the subprocesses depend on the number of gapped fermions that participate in the 
process. In Fig. E] wc have shown these rates, namely TaiA, A, A), rd(A, 0, A), TaiO, A, 0), and rd(0, 0, 0). In Fig. El 
we show four different bulk viscosities, assuming (fictitiously) that the bulk viscosity solely depends on each of these 
rates. The different suppression of the rates as a function of temperature leads to different temperatures at which the 
bulk viscosity becomes maximal (see Sec. IV Bl for a further discussion of the maxima of the viscosity). What do we 
conclude for the 2SC bulk viscosity from these curves? Or, in other words, how do these curves "add up" to yield the 
total bulk viscosity as a function of temperature? To answer this question, remember the analogy between periodically 
expanding and contracting quark matter and an electric circuit with alternating voltage, see Sec. |TT1 In that section 
we explained that one may think of a "capacitance" that behaves like the inverse of the rate 1/A, see Eq. ([T6| . If 
the rate consists of a sum of several distinct subratcs (as is the case in 2SC quark matter), small subrates can be 
neglected if there is a large, dominating subrate. The resulting "capacitance" is then determined by the dominating 
subrate, just as in an electric circuit with capacitors Ci, C2 in series where 1/Ctotai = 1/Ci + 1/C2- Therefore, the 
bulk viscosity of 2SC quark matter is dominated by the subprocess of ungapped fermions, Td{0, 0, 0), although each of 
the subprocesses with gapped fermions would, for certain temperatures and in the absence of Td{0, 0, 0), yield a much 
larger bulk viscosity (as can be seen by comparing Fig. [7] with the final results. Fig. [8]). Note that it is crucial for this 
argument that all subprocesses (all "capacitors") affect one single chemical nonequilibrium (one single "current"), 
characterized hy = fig — f-id- 

In Appendix IB] wc show how this argument is altered when there is more than one independent S^. To illustrate 
this point let us assume for a moment that the strong interactions, which can change the color of the quarks, are 
no faster than the weak interactions. In this case, we must not consider the "accumulated" process u + d <-> u + s, 
but nine processes iu + jd <-*■ ju + is, where i and j each stand for one of three colors, i, j = red, green, blue. And, 
importantly, there is not only a single Sfi, but nine (5/x„'s, one for each of these subprocesses. In this picture, the 
bulk viscosity is sensitive to the re-equilibration of each of the nine out-of-equilibrium quantities (5/x„, leading to a 
significantly different result. Roughly speaking, the result would be a curve that follows all the maxima of the curves 
in Fig. [71 However, this picture is unphysical since it ignores the effect of the strong interactions. They provide 
color-changing processes on a time scale much faster than the weak processes. We show in Appendix [B] that, after 
taking into account the strong processes, one arrives at the simpler picture that has been used in Sec. |TT1 There is 
only a single out-of-equilibrium quantity 6fi and the subrates add up like capacitors in series in an electric circuit. 
Then, the physical result for the total bulk viscosity is given by Fig. [51 which we discuss in the following. 
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FIG. 7; (Color online) Bulk viscosities calculated as if the only equilibration processes were those that contribute to the 
rd(0, 0,0) subrate (black), the rd(0, A,0) subrate (blue), the rd(A,0, A) subrate (red), and rd(A, A, A) subrate (green). We 
fixed uj/{2tv) — 10^s~^ and jj, — 500 MeV. The temperature scale is linear. 



B. Bulk viscosity of 2SC quark matter 

In Fig. [8] we show the bulk viscosity for unpaired and 2SC quark matter for two different oscillation frequencies 
Lo /{2ti) ~ 200 s~^ (left panel) and lo /{2ti) = 1000 s^^ (right panel). We now explain the features of these plots. 
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• The bulk viscosity peaks at a particular temperature. This is expected from Eq. ([21]) : the bulk viscosity as a 
function of the rate 7, for a fixed external frequency w, has a maximum at 7 = tj, when the frequency of the 
externally imposed compression matches the microscopic equilibration rate. Since the microscopic rate ^{T) is a 
monotonic function of the temperature (see results of the previous section), we see this maximum when we compute 
the bulk viscosity as a function of T. 

• The peak viscosities are easily determined by setting 7 = in Eq. (j2ip , 

" 2CZ ' ^^^^ 

i.e., they are independent of the microscopic equilibration rate. In particular, this simple relation explains the 
larger maximum values of the viscosity in the left panel of Fig. [S] compared to the right panel because of the 
different values of uj. In general, Cmax is different for unpaired and 2SC quark matter because the values of B and 
C differ from phase to phase, sec TableHl For the chosen values of to^ and A, the peak values happen to be almost 
identical. 

• The 2SC phase has its maximum viscosity at a higher temperature than the unpaired phase; near and above this 
temperature (up to T ~ 2Tc/3), the 2SC phase has a larger viscosity than the unpaired phase. This follows from 
the fact that 7(T) is monotonically increasing, but the 2SC phase has a slower equilibration rate. At very low 
temperatures both rates are slow, and the higher rate is closer to w, so the unpaired phase has higher viscosity 
[4^ . At higher temperatures, both rates are above uj, but now the slower rate is closer to uj, so the 2SC phase has 
higher viscosity. 

• At low temperature, up to temperatures where the viscosity of the unpaired phase peaks, the unpaired and 2SC 
viscosities differ by a constant factor (the two lines are parallel on our logarithmic plot). This follows from C 7 
for 7 ^ and the results of the previous section where we have shown that for sufficiently small temperatures, 
roughly T < 0.3 Tc (sec Fig. [6]), the rates simply differ by a factor 9 because the contribution of the gapped modes 
can be neglected. Also for temperatures larger than the temperature at which the viscosity of the 2SC phase peaks 
up to r ~ 10 MeV, the viscosities differ by a constant factor. However, they are now in reversed order since 
^ cx 1/7 for 7^0;. 

• At the critical temperature, here chosen to be Tc = 30 MeV, the rates of the two phases become the same (the 
2SC gap vanishes). However, due to the "locking" of Fermi surfaces in the 2SC phase and the resulting difference 
in the coefficients entering the viscosity, the viscosities are different. From the values in Table H] we conclude 
Cunp(T'c)/C2Sc(T'c) = 9/4 + 0(m^//i^), in agreement with the results shown in Fig. [51 
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FIG. 8: Bulk viscosities as a function of temperature for 2SC phase (solid line) and unpaired phase (dashed line), for a chemical 
potential /i = 500 MeV, and oscillation frequencies uj/{2-k) — 0.2 ms~^ (left panel) and uj/{2n) = 1.0 ms~^ (right panel). Note 
that temperature is plotted on a logarithmic scale. For 30keV < T < 20 MeV and 70keV < T < 20 MeV, respectively, 2SC 
has a higher bulk viscosity. 



VI. CONCLUSIONS 

We have computed the bulk viscosity of three- flavor color-superconducting quark matter in the 2SC phase originating 
from the nonlcptonic weak process u + d u + s. To this end, we have presented a detailed calculation of the rate for 
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this process for arbitrary temperatures. Since 2SC quark matter contains ungapped strange quarks as well as gapped 
and ungapped up and down quarks, this process contains several subprocesses. Each of these subprocesses that 
contains at least one gapped quasiparticle mode can be further divided into "subsubprocesses" which are similar to 
the ones in ordinary superfluids such as '^He. In both color superconductors and ^He, these subsubprocesses allow for 
creation and annihilation of quasiparticles from the condensate. Therefore, the process u + d <-*■ u + s has contributions 
from decay and coalescence processes. We have presented a formal derivation of the total rate, starting from a kinetic 
equation and making use of the microscopic fermion propagators, resulting in an expression that naturally contains 
all these contributions. Our results arc not only interesting for the special case of the 2SC phase but also, because of 
the multitude of subprocesses, provide a general quantitative analysis for the effect of a fermion condensate on the 
considered nonleptonic process. 

We have evaluated the rate numerically, observing the expected exponential suppressions for small temperatures. 
For subprocesses with three gapped modes we find a suppression of the form e"^'^/"^, while two or one gapped mode(s) 
yield e"'^!^ . Upon evaluating the collision integral for temperatures up to the superconducting transition temperature, 
we have also studied the temperature regime where the exponential functions are no longer valid approximations for 
the rate. 

The results for the rate, which translate into a typical time scale of the weak process, have been used to compute 
the bulk viscosity as a function of temperature for a given external oscillation frequency, chosen to be of order 1 kHz, 
which corresponds to the fastest rotation rate observed for pulsars. We have pointed out that the bulk viscosity is the 
response of the quark system to this external oscillation just as an electric circuit responds to an external alternating 
voltage by an induced current. In this sense, we have identified a "resistance" and a "capacitance" of the quark 
system in terms of the microscopic process. Using this analogy, it is easy to understand that the behavior of the bulk 
viscosity assumes a maximum when the microscopic time scale matches that of the external frequency. 

Our calculations show that, for a transition temperature of = 30 MeV and oscillation frequencies in the kHz 
range, the bulk viscosity of 2SC quark matter is larger than that of unpaired quark matter over a wide temperature 
range. This range is approximately given by Vd~^Tc <T < 2Tc/3 for an oscillation period of one millisecond. For 
longer oscillation periods, this range is even larger. As discussed above, this is because in that temperature range 
the slower equilibration rate of 2SC matter resonates better with kHz frequency oscillations than does the faster 
equilibration rate of unpaired quark matter. This result has potential physical implications for the behavior of the 
star early in its life, before its internal temperature drops below about 0.1 MeV. It is also interesting to note that at 
these early times neutrinos are trapped, and the resultant Icpton number chemical potential can favor a 2SC phase 
with strange quarks [ll| (but see also [il]). Moreover, bulk viscosity is likely to be the dominant factor for damping 
of potentially unstable r-modes in this very temperature regime [2^ [23]. Although the star cools down below 0.1 
MeV within days, this short time might potentially be long enough for instabilities to grow and spin down the star 
drastically [4^ . For hot newly born stars, we have thus found that quark Cooper pairing actually works in favor of 
damping of unstable modes. For much smaller temperatures, T ^ 0.1 keV, shear viscosity effects become important, 
and Cooper pairing is expected to reduce the shear viscosity and thus the damping of the unstable modes [431 ■ We 
have found a similar effect for the bulk viscosity: for sufficiently small temperatures, T < 10~'^Tc, the bulk viscosity of 
2SC quark matter is smaller than that of unpaired quark matter. Besides the damping of unstable rotational modes, 
bulk viscosity is also an important quantity for time scales of radial pulsation of compact stars [45j . Our results can 
serve as an ingredient for calculations addressing the viscous damping of these pulsations. 

Finally, we emphasize that the results for the bulk viscosity in the 2SC phase cannot be extrapolated to the 
CFL phase. In this phase, all fermions acquire a gap in their energy spectrum and thus the rate of the process 
M + c? u + s will be smaller than in the 2SC phase. In particular, this rate will be exponentially suppressed for small 
temperatures. However, kaons (which have nonzero strangeness) serve as a more effective source for re-equilibration, 
and their contribution to the bulk viscosity has to be taken into account (46j . 
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APPENDIX A: BULK VISCOSITY INCLUDING LEPTONIC PROCESSES 

In this appendix, we derive an expression for the bulk viscosity which is valid for a system that not only converts 
s into u quarks via the non-leptonic process considered in the main part of the paper, but also contains leptonic 
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u + d 


u + s , 






u + e 


d+v , 


r2 




u + e 


s + V , 







processes. We shall consider the following three processes, employing the linear approximation for the several rates 
as in Sec. HIl 

(Ala) 
(Alb) 
(Ale) 

The quantities Fi, Ai, and correspond to Td, A, and (5/x in Sec. [Ill In chemical equilibrium, i5/ii = 6^2 — — 0. 
Out of equilibrium, the quantities 

6fii = /.is - fid , (A2a) 

6fJ-2 = fid - /i« - yUe , (A2b) 
Sflz = fls - flu - fJ-e= Sfll + 5fl2 , (A2c) 

are nonzero (we assume fii, = Q). Note that there are only two independent Sfi's. According to the above microscopic 
processes, the densities for electrons and quarks change as follows. 



T2 + = A3(5/ii + (A2 + Xz)5fi2 , (A3a) 
Fi - r2 = Xi5fii - >^25fi2 , (A3b) 



dUe 

~dt 

dud 

~dr 

dus dud due , , 

^ = -^'-^' = -^'^^ 
dn„, dn„ 

We can thus express all density changes in terms of the changes in rig and nd- Analogous to Sec. |TT1 we write the 
pressure as 

- p' + ( f ) + (If) + (If) + (If) + (If) 

= P0+ (^) SV + {Bi + B2)5n, + BiSna , (A4) 



where 



Si ^ ( nd^J^\ ) , (A5a) 



dud / V 



- ( nu'^) + f..^) - fn.|^) , (A5b) 



and we have made use of the relations between the various densities given in Eqs. (jA3[) . For the purpose of this 
appendix, it is sufficient to consider unpaired quark matter. In the 2SC phase, Bi and B2 would have additional 
terms due to the pairing of u and d quarks. Since the conclusions of this appendix do not depend on the values of 
Bi and B2 we may for simplicity proceed without taking into account the effects of pairing. Then, the average power 
dissipation, cf. Eq. ([§]), is 



ldW\ 

ydT/ 



ai {5fii{t)5v{t)) + 02 {5fi2{t)5v{t)) , (A6) 



where we abbreviated 



ai = BiXi + {Bi + B2)X3 , (A7a) 
02 = B2A2 + (Si + B2)A3 . (A7b) 



The changes in the chemical potentials are derived as in SeclTTJ cf. Eqs. (fTT|) and (|12p . 



—g^ = Bi-^-ai5fii-f3i6fi2, z = l,2, (A8) 



20 



where 



ai = Cj Ai + C,f A3 , 

A = C,^(A2 + A3) - Cf A2 , 



with 



\dndj f 



C^2 = 



duu 



drip 



(A9a) 
(A9b) 



(AlOa) 
(AlOb) 



The single differential equation thus has become a system of two coupled differential equations (|A8p . With the 
ansatz = Re(5/x^e"^*) wc may reduce the system of differential equations to an algebraic set of four linear equations 
for the four variables Re 6^1 , Im Sfj-i , Re (5^2 > I^^^ ^^^2 j 



ojUcSji^ = LuBi Svo ^ oiilmSfii ~ (3ilmSiJ.2 , 
wlmJ/^" = tti Re (5/x° + /3i Re (5/i2 • 



(Alia) 
(Allb) 



From Eq. (jA6p and Sv{t) = Svocosut we conclude that only the real parts of the complex amplitudes (5/U° enter the 
power dissipation. They are given by 



Red Hi ~ 6vquj 



Si(/3| 



a2/3i) - S2/3i(ai + /32) + -Biw^ 



(a2/3i - ai/32)2 + (2a2/3i + af + /3|)tj2 + ^4 



6vo oj 



B2{a\ + a2Pi) - Bia2{ai + (32 



B2L0^ 



Consequently, Eq. 



(a2/3i - ai/32)2 + (2a2/3i + a? + /3|)u;2 _^ ^4 
6| and the definition ^ yield the bulk viscosity 



C = ai 



Re Re,5Ai§ 



02 



P2C.' 



(Al2a) 
(A12b) 

(A13) 



where 



P2 

Qi 
Q2 



(A1A2 + A1A3 + A2A3) [bI {{CI - CI)[ClX2 - C2^(A2 + A3)] + [CI - )(C|Ai + C2^A3)} 

+ B1B2 {{C^ - CI - C|)[Cf Ai + (C2" - C2^)A2] + [Ct{Cf - CI + C|) - + CI + Ci)]X^] 

+ Bl[{CifM + (Ci^- - Ci)Ci\2 + (Cf + C|)Ci-A3]) , (Al4a) 

BlXi + b2A2 + (Bi + S2)2A3 , (A14b) 
(A1A2 + A1A3 + \2\:^?{ClCi - CiC'^f , (Al4c) 
(Cf Ai + Ci^A3)' + 2(C|Ai + C|A3)[Ci^(A2 + A3) - A2] + [C|A2 - C2^(A2 + A3)]' . (A14d) 

With the densities (of unpaired quark matter) n„/d — A^^/^/tt', lis — {^J■l — ml)^^'^/Tr^, and Ue = fil/iSn'^) we have 

(A15a) 
(A15b) 
(Al5c) 

From the definitions it it clear that Bi and Cf correspond to B and C in SccUll respectively. The result (jAlSp together 
with the definitions (|Al4p shows that the three processes are entangled in a complicated way. Let us consider some 
limit cases in which the processes disentangle. In particular, we are interested in the case Ai 3> A2, A3, since we expect 
the rates for the leptonic processes to be much smaller than the rate for the non-leptonic process. In this case. 




Xi{X2 + X3){Cf - Ct)Cf + Lo^ 



A?(A2 + A3)2(CfC| - CfCiy + (Cf Ai)2c^2 + ^ 



(A16) 
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order of time scales 


c 


t<J > 71 2> 72 


7i/cj2 


Lj ~ 71 > 72 


71/(7? +0.2) 


71 S> S> 72 


72/t<J^ + 1/71 


71 » li) ~ 72 


72/(72 +'^^) 


71 3> 72 > (X) 


1/72 



TABLE III: Qualitative behavior of the bulk viscosity for different limit cases in which the time scales given by the external 
oscillation frequency u, the non-leptonic process 71, and the leptonic process 72 are clearly separated. 



In order to obtain some qualitative insight we have to compare the inverse time scale m with the inverse time scales set 
by Ai, A2, and A3. For simplicity, let us set A2 = A3. As a rough estimate, we shall assume that the factors containing 
one or several C's all contribute by a factor l/^l, i.e., Cf - Cf ~ Cf l/fi^, CfCl ~ CfCI - {Cf)^ Then, 
the relevant inverse time scales 

71 = ^ , 72 = ^ (A17) 

are, for fixed [is-, indeed set solely by the corresponding rates of the non-leptonic and leptonic processes, respectively, 
and we can write 

7172 

7172 + 71^^ + ^ 

Obviously, it is incorrect to neglect the much slower process in general. Which process can be neglected rather depends 
on ll). In particular, wc may consider several limit cases in which the time scales clearly separate, see Table Hill The 
general and plausible conclusion from this table is that the bulk viscosity is dominated by the process whose time 
scale is closest to the time scale set by the external oscillation. (Note that this is in contrast to subprocesses which 
contribute to one single b\i^ as considered in the main part of the paper. In this case, the fastest subprocess dominates 
the bulk viscosity even if there is a slower subprocess which, in the absence of the fast one, would produce a much 
larger bulk viscosity.) More specifically, the leptonic processes can be neglected in the high frequency limit (first row 
of Table [111)) or when the frequency is close to the rate of the non-leptonic process (second row). When the frequency 
is smaller than the rate of the nonleptonic process, the leptonic processes can only be neglected if its rate is much 
smaller than the frequency and uj/^i 2> 72/w, i.e., if the frequency is much closer to the rate of the non-leptonic 
process than to the rate of the leptonic processes. 



APPENDIX B: BULK VISCOSITY INCLUDING STRONG PROCESSES 



In this appendix, we show that the expression for the bulk viscosity derived in Sec.|TT]is only correct in the presence 
of strong processes which are infinitely fast compared to the weak process u + d u + s. The treatment in Sec. |TT] 
apparently simplifies the situation because the process u + d ^ u + s ui fact includes several combinations of quark 
colors. In principle, one has to treat all possible processes in a way analogous to the method presented in Appendix 
Rl i.e., wc have to consider the processes 



bu - 


'rbd ^ 


bu 4 


-bs, 


Ti = Ai J/ii , 


(Bla) 


gu- 


>rbd ^ 


bu 4 


- gs , 


r2 = A2 6^2 , 


(Bib) 


ru - 


'rbd ^ 


bu 4 


- rs , 


r3 = A3 J/i3 , 


(Blc) 


bu H 


-gd ^ 


gu- 


\-bs, 


r4 = A4 (5/^4 , 


(Bid) 


bu - 


- rd ^ 


ru - 


hfes, 


^ A5 6fi5 , 


(Blc) 


gu^ 


- gd ^ 


gu - 




Te = Ag , 


(Blf) 


ru - 


- gd ^ 


gu - 


1- rs , 


Ty = A7 S^l7 , 


(Big) 


gu - 


- rd ^ 


ru - 


^ gs , 


Fg = Ag 6 fig , 


(Blh) 


ru - 


- rd 


ru - 


h rs , 


Fg = Ag SfiQ . 


(Bli) 



As in Eqs. (jAip we assign a rate F^ to each of the processes and assume these rates to be proportional to the 
corresponding J/i.;. There are nine different processes, because we can assign one of three colors (r, g, b) to each of the 



22 



two (color-conserving) vertices of the process u + d ^ u + s. From the results in Secs. lIIII and llVl we know that not all 
of the subprocesses have different rates. Due to the remaining SU (2) color symmetry, we rather have four different 
rates depending on the number of gapped quasiparticles participating in the process. Using the notation introduced 
in Eq. ([55)) . we have 

Ti = rrf(0,0,0), (B2a) 
T2 = = rrf(0,A,0), (B2b) 
r4 = Tg = rrf(A,0,A), (B2c) 
Te = = Tg = Tg = TdiA, A, A) . (B2d) 

However, for the purpose of this appendix we do not have to make any assumptions for the rates and may treat all 
nine rates, given by Ai, . . . , Ag, as independent. 

The various 6fii^s can be written in terms of nine chemical potentials (one for each combination of color and flavor). 
It turns out that only five of them are independent, 

(B3a) 
(B3b) 
(B3c) 
(B3d) 

(5/i7 + S^g , (B3e) 

(B3f) 
(B3g) 
(B3h) 
(B3i) 

We may thus choose 5^2, ^Me, S^7^ ^1^9, as independent variables. The density changes due to the above processes 
can be read off from Eqs. (|Bip . 





= fJ-bs — 


^J■bd , 
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fJ-gs — 
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— I^^bd , 
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f^rs 




— l^bd = 


S^J-2 


- S^ie H 


-Sfi7, 




f^gu ^ 


- M6s - 


M6u 


- fJ.gd = 




- (5/i2 - 






Mru 


- M6.S — 


M6u 


— fJ-rd = 


S^ii 


- Sfi2 H 


' Sfl6 - 




= fJ-gs — 
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f^gu ^ 


fJ>rs 
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— f^ru 


- fJ-gs - 




— fJ-rd — 




- (5/i7 - 


f dllg , 


(5/ig 


— /^rs 


fJ-rd ■ 













^ ^ Ts-T, + T7~Ts, ^^T. + Ts+Tg, ^ ^ -Pg - - Tg , (B4a) 

at at at 

^ ^ T2-T,-r7 + Ts, ^^r, + r, + T7, ^^-r^-Te-rg, (B4b) 

at at at 

—rr = - + r4 + , — — = Ti + r2 + Ta , — — = -Ti -Vi-T^. (B4c) 

at at at 

Obviously, the total change of the u flavor density vanishes, while the total change of d and s flavor densities has the 
same absolute value but opposite signs, 

duu dnd dus ^ , . 

1=1 



where 



' n^i , i = u,d, s . (B6) 



If the bulk viscosity was determined solely by the above nine subprocesses of the weak process u + d^u + s,we would 
have to set up a set of five coupled differential equations for the variables S^i, 5^2, <^M6, Sfig, analogous to Eq. 
(|A8p . Its solution would yield the bulk viscosity through an equation analogous to Eq. (|A6|) . We have checked that this 
procedure would yield a result that differs from the one presented in Sec. |Vl In particular: While the bulk viscosity 
given by Eq. ([2T|) is dominated by the fastest of the subprocesses (here Ai, . . . , Ag), the setup presented above would 
potentially allow all subprocesses to contribute significantly to the bulk viscosity (the dominant subprocesses would be 
determined by the external oscillation frequency to). In other words, this setup requires all chemical nonequilibria to 
re-equilibrate separately, and the one that resonates best with the external oscillation gives the largest contribution. 
In the remainder of this appendix, however, we show that this picture is invalidated by the presence of the strong 
interactions which allows the fermions to change their color on a much faster time scale than the one given by the 
weak interaction. 
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The relevant strong processes change the colors of two different quark flavors. Since there are three pairs of distinct 
flavors and three pairs of distinct colors, we have to consider the nine processes 

ru + gd gu + rd , ru + gs gu + rs , rd + gs ^ gd + rs , (B7a) 
ru + bd hu -\- rd , ru + bs bu + rs , rd + bs bd + rs , (B7b) 
gu + bd ^ bu + gd , gu + 6s bu + gs , gd + bs ^ bd + gs . (B7c) 

Assuming these processes to be infinitely fast (this is equivalent to the assumption that the system be in chemical 
equilibrium with respect to all of these processes at any time) translates into the following conditions for the chemical 
potentials, 

l^ru ^gd l^gu P'rd 5 l^ru ~^ ^gs ^^gu l-^rs 7 f^rd ~^ ^gs f^gd /^rs 7 (BSa) 

+ Mfcd ~ Mbu ^ Mrd = , /i„i + Hbs — M6u ^ Mrs = , fird + Mbs — Mbd — Mrs = , (B8b) 

Mgii + Mfcd ^ Mbu ^ Mffd = I Mffu + Mfcs ^ Mbu ^ Mffs = , + fibs — Mbd ^ Mgs = . (B8c) 

These equations result in four independent constraints. Consequently, we may express four chemical potentials in 
terms of the other five, e.g., 

^gd Mrd ^ f^gu Mrit 7 Mps Mrs f^gu Mru ; (B9a) 

Mfcd = Mrd + M6m ^ MrM , Mfcs = Mrs + Mbn ^ Mr« ■ (B9b) 



Inserting these equations into Eqs. (|B3[) reduces the mmrber of independent (5/i's to one, 

(5/^1 = ... = (5^9 = /i^s - ^rrf = . (BIO) 

In order to set up the differential equation for this 5fi, we have to know the dependence of the densities on the remaining 
chemical potentials (after eliminating /ig^, fj.bd, fJ-gs, fJ-bs via Eqs. (|B9[) ). Then, in the 2SC phase, the densities of the 
paired quarks ti^u, nrd^ rigu, rigd depend on the average chemical potential {^rd + tJ-gu)/'2- The densities of the unpaired 
quarks are straightforwardly obtained with Eqs. (|B9p . Let us, for simplicity, assume that fi^u = Mgw = Mfcu = Mm- 
Then, we use the fiavor densities from Eq. (jB6P to obtain 



d6iJ, f dSfi\ dSv f d5fj,\ dn. 



l — U 



with drii/dt given in Eq. (|B5p . The same procedure is applied to the expression for the power dissipation, cf. Eq. 
([9]). Upon using li + . . . + Ig = A with A defined in Eq. ([8]), we see that the final result is identical to Eq. (|2ip . 
Consequently, we have proven that, due to (infinitely fast) strong processes, one can ignore the subtleties related to 
different J/x's corresponding to subprocesses oi u + d ^ u + s. One rather may solely consider the rate which is 
the sum of the rates of these subprocesses and a single related to this rate. This is the method used in the main 
part of the paper. 



APPENDIX C: DIRAC TRACES AND CONTRACTIONS 

From the definitions (l34t and (l43l) we find 



and 



rOO(k,p) = 2(l + k.p) , (Cla) 

T*(k, p) = 2 [k' +f + tCkx pY] , (Clb) 

T'°(k,p) = 2[k'+f~i{kxpy], (Clc) 

r'J'(k,p) 2[5'i{l-k-p) + k'p^ +yf + ie'^^{k'^^p^)], (Cld) 

Z^°°(k,p) = -2(l + k-p) , (C2a) 

(k, p) = 2 [k' +f + iikx pY] , (C2b) 

Z^*(k,p) = -2[k'+f-i{kxpY], (C2c) 

U'^ (k, p) = 2 - k • p) + k'f + + i e'^'^ {k^ - /)] . (C2d) 
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Consequently, we obtain the contractions 

T''''(P4,P1)V(P3,P2) = 16(l-p4-P3)(l-Pi •P2), (C3a) 
^''''(Pl, P4) V(P3, P2) = 8 {(1 ~ P4 • P3)(l - Pi • P2) + (1 + Pl • P3)(l + P4 • P2) - (1 + Pi • P4)(l + P3 ' P2) 

- i [(P3 + P2) ■ (Pl X P4) + (pi + P4) • (p2 X P3)]} , (C3b) 
W^''(P4, Pl) V(P3, P2) = 8 {(1 - P4 • P3)(l - Pl • P2) + (1 + Pl • P3)(l + P4 • P2) - (1 + Pl • P4)(l + P3 ' P2) 

+ i [(P3 + P2) ■ (Pl X P4) + (pi + P4) • (p2 X P3)]} . (C3c) 

The imaginary terms in Eqs. (jC3b[) and (jC3cP do not appear in the colhsion integral since they vanish after angular 
integration, see Appendix ID 21 

APPENDIX D: ANGULAR INTEGRATION 

In the first part of this appendix we compute the angular integrals in Eqs. and (|49p . In the second part we 
show that the angular integral over the imaginary terms in Eqs. (|C3b[) and (jC3cp vanishes. 

1. Angular integration in Eqs. (|48|) and (|49p 

We have to compute 

A' = y" d^i j (1^2 J dfls J dn^l - Pl • p2)(l - P3 • P4)'5(pi + P2 - P3 - P4) ■ (Dl) 
One may start with the d^l2 integral and rewrite the (5-function in terms of a modulus and an angular part 

S[P2 - (P3 + P4 - Pl)] = \S{P2 - |P3 + P4 - Pl|)'5(f^2 " ^) , (D2) 

P2 

where il stands for the angles of the vector ps + p4 — pi. Thus, the effect of the angular 5-function is to replace the 
direction p2 with (p3 + P4 — Pi)/|p3 + P4 — Pi| (the numerator being actually P2 because of the radial (5-function). 
Consequently, with P = p3 + p4 we obtain 



K = \ I dfli I dQs I d^i 



l--(pi-P-pi) 

P2 



(1-P3-P4)%2-|P-P1|). (D3) 



Next, we may straightforwardly perform the dQi integration by choosing the z-axis to be parallel to P. As in Ref. [40| 
we denote 

P^J = \Pi - Pj h Ptj = Pl + P] ■ (D4) 

Then, 



dVLi 



l--(pi-P-pi) 

P-2 



5{p2 - |P - pil) = ^^Q{P - Pi2)e(Pi2 - P){Pl2 - P^) ■ (D5) 



PiP 

To obtain this result, note that 8(a — \b\) ~ Q{a — h)Q{a + 6). The remaining angular dependence is reduced to the 
angle between p3 and p4, present in P. Hence, one integration, say d^^, is trivial while we are left with the azimuthal 
integral from (if23. 



K =^ 
P1P2 



/•I 1 

j ^ dx -e(p - pi2)e(Pi2 - P){pf2 - ^')(i - 

/ dPQ{P-p,2)Q{Pl2 - P){P!2 - P^){Pk - P' 

Jp34, 



pIpIpIpI Jp, 



(D6) 



Since Ptj > ptj, there are four orders of pi2,P34, -P12, P34 that yield nonvanishing integrals, 

P12 < Pm < P12 < P-ii , (D7a) 

P3i < P12 < P12 < Psi , (D7b) 

P34 < P12 < P34 < P12 , (D7c) 

P12 < P34 < P34 < P12 ■ (D7d) 
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The two other possibiUties, pi2 < P12 < Psi < -P34, and ^34 < P34 < pi2 < P12 lead to vanishing integrals. 
Consequently, 

^ = 2222 HP12, A2,P34, ^34) , (D8) 

where we defined 

L(a, b, c, d) = e(c - a)e{d - b)e{b - c) J{c, b, b, d) + Q{a - c)Q{d - b) J{a, b, 6, d) 

+ e(a - c)e(6 - d)e{d - a) J{a, d, b, d) + Q{c - 0)6(6 - d) J{c, d, b, d) , {D9) 

with 

i-b 



J{a,b,c,d) = f dP [c^ - P^){d'^ - P^) 

J a 



c^d^ib - a) - -(c^ + d^){b^ ~ a^) + 1(6^ - a^) . (DIO) 
3 5 



With a simple relabeling of momenta we obtain the result for the other two terms in Eq. (|49|) 

47r3 

piphlpz 

and 

47r3 
pfp^PlPt 
Consequently, with the definitions 

I{pi,P2,P3,Pa) = L{pi2,Pi2,Pu,Pu) , (D13a) 

Hpi,P2,P3,Pa) = L{pi2,Pi2,P34,P34)+L{p24,P24,Pl3,Pl3)+L{pi4„Pl4,P23,P23), (D13b) 

we arrive at Eqs. ((50|) and ([5T|) . 



[ dfli [ dn2 [ dfl-i [ dVti{l + pi ■ P3)(l + P2 ■ P4)^(Pl + P2 - Pa - P4) = 2 ^2^ 2 2 -^(P24, f24,Pl3, Aa) , (Dll) 

J J J J PlPmPi 
id 

[ dfli [ dn2 [ dn^ [ dni{l+Pi-pi){l+P2-p3)5{Pl+P2-P3-Pi) ^ ^^^-2-^L{PU,PlA,P23,P23)- (D12) 
J J J J P1P2P3P4, 



2. Angular integration over imaginary terms in Eqs. (|C3bp and (|C3cp 

Here, we show that the integral 

/o = y" dill J d02 J dils J dn4[{p2 + Ps) ■ (Pi x P4) + (pi + P4) • (p2 x P3)](5(pi + P2 - P3 - P4) (D14) 

is zero. After integrating over one of the angles, say over O2, the scalar products containing a cross product can all 
be expressed in terms of one single scalar product. 



P2 



Pi + P3 + Pi 
P2 



j dfli j dn^ y"(if24Pi • (p3 X P4)(5(p2 - |P3 + P4 - Pil) ■ (D15) 



Let's consider the next integration, say over ili. The (5-function ensures that the distance of pi to a fixed vector, 
namely P3 + P4, equals a fixed value, namely p2- So all allowed vectors pi sit on a sphere around P3 + P4 with radius 
P2- The value of the integrand is given by the projection of pi onto another fixed vector, namely p3 x p4. Since this 
vector is orthogonal to ps + P4, the integral over all these projections vanishes. 

For the formal proof, one uses a frame with z-axis parallel to P3 x p4 and an x-axis parallel to ps + p4 (as above, 

P = P3 +P4), 



/ 



dilipi ■ (p3 X P4.)5{p2 - |P3 + P4 - Pl| 



dO sin 6* cos 6* 
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dip5{p2 



P"^ + p\ — 2Ppi sin 9 cos (p) 



2P2 

Ppi Jo 



dO- cos 

I sinwl 



pI + p^-pV ' 

2Ppi sine 



-1/2 



0, (D16) 



because the integrand is odd in 6 over the integration range. 
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<5^ < Ao 


Ao <5^i < 2Ao 


2Ao < (5/1 < 3Ao 


3Ao < Sfi 


rd(A,A,A) 


g-(2A-iM)/T 


e-A/T 


g-(3A-4M)/T 


1 


rd(A,0,A) 


e-A/T 


g-(2A-«M)/T 


1 


1 


rd(o,A,o) 


g-{A-iM)/T 


1 


1 


1 


rd(o,o,o) 


1 


1 


1 


1 



TABLE IV: Leading small temperature behavior of the rates with three, two, one, and zero gapped modes. A "1" indicates 
that the rate does not vanish at T = 0. 



APPENDIX E: RATE FOR SMALL TEMPERATURES 



In this appendix, we determine the leading contribution for the rates Td{0, A, 0), rd(A, 0, A), rd(A, A, A) at small 
temperatures. To this end, we analyse the following expression which originates from the Fermi distributions in the 
integrands of Eqs. ([50]) and (|5ip (prcfactors arc ignored since we only determine the leading contribution), 

rd(Ai,A2,A4) - J2 (e'=iV^?+'^? + l)"^ (^^e,^-4+vi + 1^ 



616264 
X 



where Xi = {pi — ^i)/T, ipi = A(T)/T, y = 6fi/T. This expression is integrated over xi, X2, X4, each ranging from 
—00 to 00 (the lower boundary —oo actually being an approximation for —n/T). We perform the sum and, for small 
T (hence large (fi), approximate 



rd(Ai,A2,A4) 



1 1 



1 1 



+ 



^y/ xI+^pI + ^J Xl + if'l-y _|_ ^^Jxl+ifl p— y ^ ^^Jxl+ipl + ^Jxl+^pl + ^Jx'l + ipl 
1 1 



x'l+ipl + ^J xl+ifl + ^J xl+ipl-y ^ ^ ^y/xl+fl-y _|_ gy/xl+ipl+y/ xl+<fil 
1^1 
^y/^l+vl-V ^ gy/xj+'^j+y/xl+ifl ^y/xl+v^^^ + y/x^+ifil-y _j_ ^y/xl+ipl 



(E2) 



At this point, we disti nguish be tween the different numbers of gapped modes. First, let us discuss the case (fi = (p2 = 
(fis = (fi. We expand y^xf+ip^ :^ ip + xf/ {2ip) and introduce the new (integration) variables Zi = Xi/y/2(p, 

rw(A, A, A) ~ 5 5 — ^ + z — z T o — o — T ~t~ o — 9 — ^ . (E3) 

Now note that for sufficiently small T, both if and y have the same temperature dependence. 



Therefore, the limiting behavior of rrf(A, A, A) is different for the four cases < Aq, Aq < 6^ < 2Aq, 2Aq < Sfi < 
3Ao, and 3Ao < Sfi. In particular, the latter case leads to a nonzero value of F^ for T 0. This can be seen from 
Eq. (jESp . where the 1 in the denominator of the fourth term on the right-hand side survives for BAq < S^. Wc show 
the results for all four cases in the first row of Table IIVI 

Next, we consider the case ipi ~ (p4 = ip, (p2 = 0. The analogue to Eq. (jE3P is 



rd(A,0,A) 1.1.1 



g(p-yg2| _|_ ^ipgZj+X2 g2ip-yg2^+2| _|_ g93-J/gX2+2^ _j_ giPg2f 



g-y _|., g2tpg2f +a:2+2^ g2ip-i/g2j+a:2-|-2| _j_ ^ ^-y^X2 g2i^g2j+2| 



(E5) 



Again, the results for different values of 5^ are shown in Table HV 
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Finally, for ipi = ip4 = 0, ip2 = (p we have 

r4o,A,o) ^ ^ ^ 



1 1 1 , ^ 

^ 5 1 5 \ 5 . fE6) 

g-J/ gVe^i+^2+^4 gi/3-yg2;i+2^+.T4 _|_ 2 g<p-yg^5 _(_ ga:i+a:4 

Here, only for 8[l < Aq, the rate goes to zero for T ^ 0. The result is shown in the third row of Table [TVl For 
completeness, we added the trivial result that for rci(0, 0, 0) the rate is nonzero at T = for all > 0, see fourth 
row in Table [TVl 
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